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SECTION  I 
INTRODUCTION 


This  report  describes  the  results  of  an  experiment  to  compare 
and  evaluate  several  approaches  to  computing  high- resolution  wavenumber 
spectra  of  seismic  noise.  The  experiment  was  conducted  prior  to  simultaneous 
ambient  noise-field  studies  at  each  array  station  constituting  the  network  of 
stations.  Frequency-wavenumber  spectra  are  extremely  useful  for  deter¬ 
mining  the  frequency  and  velocity  structure  of  the  noise  field  existent  at  an 
array  station  and  may  be  used  to  locate  and  describe  all  the  components  of 
the  ambient  noise  field.  Principal  objectives  of  the  ambient  noise  study  were 
to  determine 


The  degree  to  which  various  network  stations 
see^bpergy  from  the  same  noise  sources 


The  mechanism  and  characteristics  of  thes** 
noise  sources 


•  The  relative  proportions  of  bodywave,  surface 
wave,  and  random  noise  present  at  each  station 


An  understanding  of  these  characteristics  is  a  necessary  pre¬ 
requisite  to  the  development  and  implementation  of  network-processing  pro¬ 
cedures  for  upgrading  seismic  event  detection,  location,  and  identification. 

The  results  of  the  network  noise  analysis  are  presented  in  a 

1 

companion  report,  Advanced  Array  Research  Special  Report  No.  6.  Included 
in  that  report  are  studies  of  the  spatial  distribution  of  coherent  noise  at  five 
array  stations  in  the  northern  hemisphere.  A  principal  tool  used  in  those 
studies  is  the  frequency-wavenumber  spectral  estimation  technique  deter¬ 
mined  best  suited  for  noise-field  analysis  as  a  result  of  the  experiment 
described  in  this  report. 
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Since  the  network  array  stations  were  all  of  small  aperture, 
conventional  techniques  for  computing  wavenumber  spectra  could  not  achieve 
the  wavenumber  resolution  required  to  resolve  and  accurately  locate  multiple 
noise  sources,  especially  bodywave  noise  sources.  The  high- resolution 
approach  to  wavenumber  analysis  has  previously  been  successfully  demon¬ 
strated  for  transient  signals;  however,  it  is  also  applicable  to  coherent 
seismic  noise  studies.  Objectives  of  this  experiment  were  to  establish  the 
validity  of  the  high- resolution  wavenumber  spectral  estimate  for  wavenumber 
analysis  of  coherent  seismic  noise  and  to  determine  those  parameters  neces¬ 
sary  for  optimum  implementation  of  the  technique  for  such  analysis. 


( 


J 
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SECTION  II 


THE  EXPERIMENT 
2 

A  recently  developed  technique  for  computing  frequency- 
wavenumber  spectra  with  very  high  wavenumber  resolution  has  been  suc¬ 
cessfully  used  to  locate  event  epicenters  through  the  computation  of  wave- 
number  spectra  from  short-period  recordings  of  their  P  v/aves.  The  tech¬ 
nique,  although  equally  valid  for  computing  wavenumber  spectra  for  noise 
fields,  had  been  little  used  for  studying  the  composition  of  ambient  noise; 
there  was  reason  to  believe  that  the  parameters  for  optimum  implementation 
of  the  technique  for  transient-signal  wavenumber  location  might  not  be 
optimum  for  the  case  of  several  rather  diffuse  and  quasi-stationary  noise 
sources. 


The  reported  experiment  was  designed  to  compare  and  evaluate 
several  approaches  to  computing  high- resolution  wavenumber  spectra  of 
ambient  noise.  A  single  noise  sample  recorded  at  CPO  on  a  day  of  intense 
and  well-defined  storm  activity  (16  October  1964)  along  the  eastern  sea¬ 
board  was  used  for  all  comparisons. 

The  CPO  station  was  selected  on  the  basis  of  extens'.ve  prior 
study  of  the  ambient  noise  there;' that  study  had  shown  the  noise  to  consist 
of  a  very  characteristic  and  time-stationary  component  plus  a  component 
directly  attributable  to  /storms  on  or  near  the  North  American  continent. 

With  strong  storms  located  at  Cape  Hatteras  (tropical  storm  Isabel)  and  on 
the  Newfoundland  coast,  the  belief  was  that  the  noise  field  at  CPO  that  day 
could  be  easily  predicted  and  that  the  prediction  could  be  considered  in 
establishing  the  validity  of  the  various  techniques  being  compared. 

The  sample  chosen  was  a  6.67-min  sample  of  noise  recorded 
by  nine  short-period  vertical- component  seismometers;  Z5  was  not  usable. 

Contract  No.  AF  33  (657)- 146  ^8 
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The  low- resolution  conventional  technique  and  several  vari¬ 
ations  of  the  high- resolution  technique  were  used  to  compute  wavenumber 
spectra  for  this  noise  sample.  Figure  1  is  a  chart  of  the  various  techniques 
and  the  methods  to  estimate  the  auto-  and  crosspcwer  density  spectra  used 
as  inpu':. 


While  many  interesting  comparisons  are  possible  from  such 
a  variety  of  approaches,  those  most  pertinent  to  the  purposes  of  this  ex¬ 
periment  arc 

•  The  conventional  wavenumber  spectrum 
computed  using  a  set  of  crosspower  spectra 
derived  by  Fourier  transformation  of  corre¬ 
lation  functions  compared  with  one  computed 
using  crosspower  spectra  derived  by  direct 
Fourier  transformation  of  array  outputs 

•  Comparison  of  high- resolution  spectra  obtained 
from  the  same  two  sets  of  crosspower  spectra 

®  Comparison  of  high- resolution  spectra  computed 
with  respect  to  different  reference  sensors 

•  Comparison  of  high- resolution  spectra  obtained 
by  two  methods  of  averaging  spectra  computed 
with  respect  to  different  referent sensors 

•  Averaged  high- resolution  spectra  compared  with 
individual  high- resolution  and  conventional  spectra 

•  Comparison  of  high- resolution  spectra  computed 
using  different  signal-to-noise  ratios  (SNR)  for 
design  of  the  multichannel  filter  which  is  the  basis 
of  the  high- resolution  technique 

•  High-resolution  spectra  computed  at  frequencies 
other  than  1  Hz  from  data  equalized  on  the  basis 

of  1-Hz  calibrations  compared  with  those  computed 
from  frequency-dependent  normalized  data 
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figure  1.  Comparison  of  Wavenumber  Spectral  Estimation  Techniques 
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SECTION  III 

CONVENTIONAL  VS  HIGH-RESOLUTION  TECHNIQUES 


Conventional  and  high- resolution  methods  are  compared  to 
establish  their  suitability  for  noise  analysis.  Figure  2  (a,  b,  and  c)  displays 
wavenumber  spectra  computed  at  0.  3,  0.  5,  and  1.  0  Hz,  respectively,  from 
smoothed  correlation  and  smoothed  direct-transform  crosspower  matrices 
for  both  techniques. 

Correlations  are  computed  for  ±5  sec  over  6.67  min  of  data 
amplitude  equalized  at  1.0  Hz,  Parzen-smoothed,  and  Fourier-transformed, 
resulting  in  a  frequency  resolution  of  0.  1  Hz.  The  data  traces  also  are 
Fourier-transformed  directly  over  the  same  data  gate  and  are  Danielle- 
smoothed  to  yield  a  0.  1-Hz  frequency  resolution.  High- resolution  wave- 
number  spectra  are  computed  using  Z3  as  the  reference  sensor  and  an  SNR 
of  0.  1. 

Each  spectrum  indicates  coherent  energy  from  the  storm 
east  of  the  station  (low-pressure  center  1,  Figure  3).  At  0.3  Hz,  excellent 
agreement  in  wavenumber  location  is  seen  for  all  methods,  with  a  marked 
improvement  in  wavenumber  resolution  with  the  high- resolution  technique. 
Although  the  storm  center  is  located  almost  due  east  of  the  station,  all 
spectra  indicate  energy  from  slightly  r  outh  of  east.  Similarly,  when  using 
correlation-derived  spectra,  agreement  is  good  in  both  azimuth  and  veloc¬ 
ity  at  0.  5  Hz.  The  conventional  technique  using  smoothed  direct-transform 
data  agrees  well  with  the  correlation-derived  spectra.  The  high- resolution 
technique  using  smoothed  direct-transform  data  locates  the  peak  slightly 
farther  south,  although  there  is  also  an  indication  of  energy  from  the  east. 

The  effect  is  somewhat  more  pronounced  in  Figure  2c,  in  which  spectra  at 
1.0  Hz  are  shown.  Again,  the  conventional  wavenumber  spectra  agree 
well,  locating  the  energy  source  slightly  south  of  east. 
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The  high-resolution  spectra  tend  to  locate  the  peak  further 

south,  with  indication  of  additional  energy  from  the  northeast.  This  effect, 

probably  caused  by  the  choice  of  reference  sensor  and/or  SNR,  is  explored 
% 

further  in  subsection  D  of  Section  IV. 

Figure  4  shows  the  CPO  array  geometry  used. 

f 


Z4 


-  F="J  1  1  =» 

0  1  km 


Figure  4.  CPO  Array  Geometry 
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Direct-transform  methods  are  computationally  simpler  and 
yield  comparable  results;  therefore,  correlation-derived  crosspower 
spectral  density  matrices  are  abandoned  in  favor  of  direct-transform- 
derived  crosspower  spectral  density  matrices. 

/  The  generally  good  agreement  of  conventional  and  high- 

resolution  wavenumber  spectra  and  of  both  with  known  noise  sources  estab¬ 
lishes  the  validity  of  the  high- resolution  approach  to  wavenumber  spectra 
of  ambient  noise.  The  slight  spectral  variations  6bserveu  for  different 
high- resolution  approaches  suggest  that  more  investigation  of  the  proper 
parameters  for  high- resolution  computation  is  in  order. 
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SECTION  IV 


COMPARISON  OF  SEVERAL  TECHNIQUES  FOR 
COMPUTING  HIGH-RESOLUTION  WAVENUMBER  SPEC1  RA 

The  various  parameters  of  the  high- resolution  technique 
investigated  —  all  of  which  have  some  degree  of  effect  on  the  wavenumber 
spectral  estimate  —  are 

•  The  means  of  compensating  for  intra-array 
inequalization 

•  The  SNR  used  in  multichannel-filter  (MCF) 
design 

•  The  method  of  estimating  and  smoothing  the 
crosspower-density  spectra 

•  The  choice  of  reference  sensor  for  MCF 
design 

Investigated  also  are  two  approaches  to  the  smoothing  of  the 
wavenumber  spectra  to  average  out  the  distorting  effects  of  inequalization 
and  nonspace-stationary  components  in  the  noise  field.  The  first  approach 
averages  the  wavenumber  power  responses  of  several  MCF's,  each  designed 
with  respect  to  a  different  reference  sensor.  The  other  approach  averages 
the  wavenumber  spectra  computed  with  respect  to  several  different  reference 
sensors.  Since  each  high- resolution  wavenumber  spectrum  is  simply  the 
reciprocal  of  the  wavenumber  power  response  of  an  MCF  designed  to  whiten 
as  a  function  of  wavenumber,  these  two  approaches  amount  to  averaging 
either  before  cr  after  taking  the  power-response  reciprocal. 

A.  EQUALIZATION  VS  NORMALIZATION 

One  method  of  compensating  for  intra-array  inequalizations 
is  to  multiply  each  trace  by  a  scalar  constant  proportional  to  the  amplitude 
response  of  the  corresponding  seismometer  as  determined  by  calibration. 
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This  method  is  simple  and  easy  to  implement,  but  the  calibration  informa¬ 
tion  is  usually  only  available  at  one  frequency  and  is  only  valid  over  a  small 
frequency  range  or.  either  side  of  the  calibration  frequency;  thus,  it  does  not 
compensate  correctly  at  other  frequencies.  As  an  additional  disadvantage, 
such  inequalization  causes  as  earth  filtering  and  seismometer-to- earth 
coupling  are  not  taken  into  account. 

Conversely,  another  method  which  requires  the  normalization 
of  crosspower  matrices  does  not  require  any  a  priori  information,  takes  all 
effects  into  account,  and  is  performed  independently  at  each  frequency.  If 
the  time  traces  are  somewhat  well  amplitude-equalized,  these  two  procedures 
should  yield  approximately  equivalent  results. 

Crosspower  matrices  are  normalized  according  to  the 

equation 


where  is  the  autopower  density  of  the  i  trace  and  is  the  measured 

,  ,  ,  .th  ,  .th  ^ 

crosspower  density  between  the  1  and  j  traces. 

High- resolution  wavenumber  spectra  are  computed  at  0.  3  Hz 
from  both  normalized  and  amplitude-equalized  (at  1.  0  Hz)  crosspower  matrices 
of  smoothed  direct  transforms.  Z3  is  the  reference  sensor,  and  the  filter - 
design  SNR  is  0.1. 

Figure  5  displays  the  wavenumber  spectra  obtained  and 
indicates  that  the  methods  compare  favorably,  although  they  are  not  com¬ 
puted  at  the  equalization  frequency.  This  agreement  implies  that  amplitude 
equalization  differences  between  channels  are  small  and  that  either  process 
can  adequately  compensate  for  them.  Normalization  is  preferable,  however, 
since  the  required  scale  factors  are  readily  available  and  no  manipulation 
of  time  traces  is  required. 
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Figure  5.  Comparison  of  High-Resolution  Wavenumber  Spectra  Using 

Equalized  and  Normalized  Crosspower  Matrices  of  Smoothed 
Direct  Transforms 


B.  EFFECT  OF  VARYING  SNR 


The  SNR  used  in  multichannel  whitening-filter  design  affects 
the  dynamic  range  and,  to  some  extent,  the  wavenumber  resolution  of  the 
high- resolution  wavenumber  spectrum. 

An  evaluation  of  the  variation  in  resolution  with  SNR  is  con¬ 
ducted  to  establish  a  proper  SNR  for  the  ambient-noise  analysis.  High- 
resolution  spectra  are  computed  at  0.  5  Hz.  Normalized,  smoothed,  direct 
transforms  of  the  entire  6.67-min  data  gate  are  used.  Frequency  resolution 
is  0.  1  Hz.  The  SNR  is  varied  from  0.  001  to  0.  5,  and  Z3  is  used  as  the 
reference  sensor.  Figure  6  displays  the  spectra. 
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For  this  noise  sample,  coherent  energy  at  0.  5  Hz  from  the 
east-southeast  had  been  detected  earlier.  This  coherent  peak  is  observed 
for  signal-to-noise  ratios  of  0.  1  and  0.  5.  As  the  SNR  decreases,  the 
dominant  peak  is  observed  to  be  resolved  into  weaker  stationary  peaks 
which  are  apparently  valid;  however,  relative  comparison  of  peak  sizes 
obviously  is  not  possible.  At  an  SNR  of  0.001,  the  coherent  energy  source 
in  the  southeast  is  poorly  located  in  wavenumber  space  and  the  spectrum  is 
believed  invalid.  The  reason  for  poor  performance  when  very  small  SNR 
is  used  is  the  almost  total  lack  of  a  constraint  on  the  MCF  wavenumber 
response  at  wavenumbers  where  no  coherent  energy  is  present. 

From  previous  studies  of  transient  signals  using  high- resolution 
3 

wavenumber  spectra,  SNR's  from  0.  1  to  0.  01  are  found  suitable.  The  pre¬ 
ceding  evaluation  indicates  valid  spectra  for  ambient  noise  using  SNR's  at 
or  above  0.  1.  The  choice  of  SNR  appears  to  have  more  effect  when  using 
smoothed  direct-transform-derived  spectral  matrices.  Similar  experi¬ 
ments  using  crosspower  matrices  of  unsmoothed  direct  transforms  indicate 
a  relative  insensitivity  to  SNR  above  0.  01;  this  apparent  sensitivity  or  in¬ 
sensitivity  to  SNR  is,  to  a  large  extent,  due  to  the  effects  of  smoothing  when 
using  direct-transform  spectra. 

C.  SMOOTHED  VS  UNSMOOTHED  CROSSPOWER  MATRICES 

High- resolution  wavenumber  spectra  computed  from  both 
smoothed  and  unsmoothed  spectra  are  compared.  Figure  7  displays  wave- 
number  spectra  at  1.0  Hz  computed  using  three  different  spectral  matrices 
obtained  by  three  different  approaches  to  the  estimation  of  crosspower  spectra. 
Crosspower  spectra  are  estimated  from  Fourier  transforms  of  Parzen- 
smoothed  crosscorrelations  and  unsmoothed  and  Danielle-smoothed  direct 
Fourier  transforms  of  the  traces.  All  matrices  are  computed  over  the  same 
6. 67-min  data  gate  and  normalized.  An  SNR  of  0.  01  is  used,  and  Z3  is  the  reference 
sensor.  Smoothed  and  unsmoothed  spectra  have  frequency  resolutions  of  0 .  1  Hz 
and  0.0025  Hz,  respectively. 
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SMOOTHED  CORRELATION 


Wavenumber  spectra  computed  from  smoothed  crosspower 
matrices  are  approximately  equivalent  and  have  several  separated  peaks 
both  north  and  south  of  east.  The  spectrum  computed  from  the  unsmoothed 
direct-transform  matrix  indicates  a  single  peak  approximately  due  east, 
corresponding  to  tropical  storm  Isabel.  As  previously  observed,  at  an  SNR 
of  0.  1,  smoothed  crosspower  wavenumber  spectra  agree  well  with  those 
obtained  from  unsmoothed  crosspower. 

Wavenumber  spectra  derived  from  unsmoothed  direct  trans¬ 
forms  seem  to  represent  a  more  stable  and  more  reasonable  estimate  of 
the  coherent  noise.  Previous  work  with  transient  signals  shows  that  cross¬ 
power  smoothing  can  virtually  wipe  out  a  wavenumber  spectral  peak  due  to 
a  coherent  transient  signal  when  the  crosspower  spectra  are  smoothed  over 
a  frequency  band  that  is  wide  with  respect  to  the  rate  of  change  of  the  cross¬ 
power-phase  spectrum.  (The  rate  of  change  can  be  very  high,  even  for 
moderate-sized  arrays.  )  If  the  noise  field  is  viewed  as  consisting  of  many 
overlapping  transients,  significant  information  could  be  lost  by  smoothing. 

If  the  noise  field  is  compose^  of  time- separated  transients  in  a  random 
noise  field,  smoothing  is  necessary  to  produce  an  interpretable  spectrum 
due  to  the  random  noise;  however,  care  must  be  taken  to  avoid  losing  the 
transient  information. 

D.  EFFECT  OF  VARIATION  IN  REFERENCE  SENSOR 

Previous  studies  show  that  significant  variations  can  exist 

in  wavenumber  spectra  computed  using  different  reference  sensors  because 

of  intra-array  inequalization  and  nonspace- stationary  components  in  the 
3 

noise  field. 


An  evaluation  of  reference  -  sensor  choice  is  performed  using 
both  smoothed  and  unsmoothed  spectra.  Figures  8  and  9  display  the  effects 
of  wavenumber  spectra  at  1.0  Hz  for  smoothed  and  unsmoothed  crosspower 


matrices,  respectively.  The  matrices  are  normalized,  and  an  SNR  or  0.01  is 


used  in  the  filter  design. 
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Figure  8.  Comparison  of  High-Resolution  Wavenumber  Spectra,  as  a  y 

Function  of  Reference  Sensor,  Using  Smoothed  Correlation- 
Derived  Matrices  * 
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Figure  9.  Comparison  of  High-Resolution  Wavenumber  Spectra,  as 
Function  of  Reference  Sensor,  Using  Unsmoothed  Direct 
Transforms 
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The  smoothed  spectra  display  considerable  variation  in  spatial 
resolution  with  choice  of  reference  sensor.  The  indicated  source  to  the  east 
is  consistent  with  conventional  and  high- resolution  spectra  presented  earlier. 
The  unsmoothed  spectral  estimates  are  apparently  more  stable  and  represent 
the  best  possible  estimate  in  this  case. 

E.  AVERAGING  OF  WAVENUMBER  SPECTRA 

Interpretation  and  source  location  are  hindered  by  the  elongated 
peaks  and  the  wandering  of  peaks  sometimes  observed  in  high- resolution  spec¬ 
tra  computed  with  respect  to  a  single  reference  sensor.  These  peaks  can  be 
resolved  by  computing  several  spectra  using  different  reference  sensors  and 
then  stacking  these  spectra, 

Essentially,  this  operation  averages  out  variations  in  the  noise 
fields  as  seen  by  different  reference  sensors  and  the  effects  due  to  variations 
of  instrument  responses  and  ground  coupling.  Spatial  resolution  is  gained  by 
the  summing  technique  accompanied  by  a  reduction  in  dynamic  range  (peak 
amplitude  to  background  level).  Peaks  become  circular  rather  than  elliptical. 

Figure  10  presents  spectra  computed  at  1.0  Hz  for  two  stacking 
methods,  i.e.,  normal  and  reciprocal.  The  normal  stacking  method  consists 
of  summing  the  wavenumber  power  responses  of  multichannel  whitening  filters 
designed  with  different  reference  sensors.  The  power  responses  of  the  fil¬ 
ters,  it  should  be  noted,  are  actually  the  reciprocals  of  the  wavenumber 
spectral  estimate.  The  reciprocal  stacking  method  uses  the  same  MCF  wave- 
number  power  responses  but  averages  the  reciprocals  of  the  individual  power 
responses.  Computation  of  each  of  the  spectra  is  from  unsmoothed  direct 
transforms.  The  crosspower  matrices  are  normalized,  and  a  0.01  SNR  is 

s 

used  in  each  MCF  design. 
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Figure  10.  Comparison  of  Two  Averaging  Techniques  for  Smoothing 
High-Resolution  Wavenumber  Spectra 


Additional  discrimination  between  coherent  peaks  and  overall 
background  level  is  achieved  by  the  reciprocal  stacking  procedure  because 
there  is  a  dependence  of  peak  resolution  on  the  relation  between  reference 
sensor  and  source  azimuth.  The  appendix  presents  a  shortcut  computational 
technique  for  each  of  the  two  stacking  methods,  eliminating  the  computation 
of  individual  filters  and  filter  responses  for  each  reference  sensor.  The 
techniques,  applicable  only  to  unsmoothed  crosspower  matrices,  represent 
dramatic  reductions  in  the  total  computations  required. 


a 
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The  validity  of  seismic  noise  analysis  using  the  high- resolution 
technique  is  demonstrated.  Several  implications  arise  that  are  pertinent  to 
the  network  surveillance  problem.  Apparently,  coherent  seismic  noise  (in 
this  case,  storm-generated)  is  sufficiently  time-stationary  over  periods  of 
at  least  several  minutes  to  permit  such  analysis.  In  addition,  p  ane-wave 
assumptions  appear  to  hold  for  this  energy.  Consequently,  high- resolution 
wavenumber  analysis  should  indicate  relative  stability  of  the  background 
noise  enhancing  the  ability  to  detect  transient  signal  arrivals  from  earth¬ 
quakes  or  explosions. 


This  behavior  may  be  rather  significant  when  assessing  the 
practicality  of  such  systems  for  seismic  event  detection  and  location.  Slow 
variations  in  the  wavenumber  spectra  caused  by  movement  or  changes  in 
intensity  of  storms  might  be  tracked  adaptively  and  those  noise  components 

H 

removed  from  the  spectrum. 


Several  conclusions  are  drawn  regarding  various  approaches 
to  computation  of  the  high- resolution  wavenumber  spectra.  Wavenumber 
spectral  estimates  made  from  unsmoothed  crosspower  spectra  appear  more 
stable,  less  sensitive  to  filter  design  SNR,  and  computationally  much 
simpler  than  those  made  from  smoothed  crosspower  spectra. 


Summation  of  several  MCF  wavenumber  power  response,  where 
each  MCF  is  designed  with  respect  to  a  different  reference  sensor,  stabilises 
the  wavenumber  estimate. 


«» 

Summing  of  individual  spectral  estimates  (i.  e.  ,  the  reciprocals 
of  the  MCF  wavenumber  responses)  yields  better  resolution  of  coherent  peaks 
relative  to  the  background  level.  0 
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APPENDIX 


DERIVATION  OF  FORMULAS  FOR  DIRECT  COMPUTATION 
OF  SMOOTHED  HIGH-RESOLUTION  WAVENUMBER  SPECTRA 

The  high- resolution  wavenumber  spectra  smoothing  techniques 
compared  in  this  report  would  both  be  prohibitively  complex  and  expensive  if 
each  individual  MCF  had  to  be  designed  and  its  wavenumber  power  response 
computed.  Fortunately,  when  unsmoothed  direct-transform-derived  cross¬ 
power  spectra  are  used,  these  smoothed  wavenumber  spectra  can  be  obtained 
directly  at  a  fraction  of  the  computational  effort. 

The  derivation  of  the  direct-solution  formulas  for  both  smoothing 
techniques  uses  a  theorem  for  the  exact  inverse  of  a  special  type  of  matrix. 

This  special  type  of  matrix  is  the  same  form  as  the  matrix  which  must  be  in¬ 
verted  in  obtaining  the  multichanr.  *1  v/hitening  filtern.  The  theorem  states 
that,  for  any  column  vectors  B  ,vi  C  and  nonsingular  matrix  A  for  which 

rn 

A  +  BC  1  is  nonsingular, 

(a  +  bSt)-1  =  a-1  -  £~1asTf'1 

1  +  £T  A_  B 

where  superscript  T  denotes  transpose. 

A.  AVERAGE  OF  WAVENUMBER  POWER  RESPONSES  OF  N  MULTICHANNEL 
WHITENING  FILTERS 

Let  X  be  an  Nxl  complex  column  vector  with  elements  com¬ 
posed  of  the  N  complex  numbers  obtained  by  direct  Fourier  transformation 
of  the  outputs  of  an  N  element  array.  The  crosspower  spectral  matrix  is 

PI 

then  XX  where  superscript  H  denotes  conjugate  transpose.  Sought  is  a  mul¬ 
tichannel  filter  which  will  whiten,  as  a  function  of  wavenumber,  the  space- 
time  field,  the  crosspower  matrix  of  which  is  XXH-  The  equation  to  be  solved 


A-  1 
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for  the  complex  Nxl  filter-weight  vector  is  given  by 


f*  =  (si  +  xxH)_1  r.  =  -  z  r. 

1  \  /is  1 


where  the  asterisk  denotes  conjugate,  s  is  a  real  scalar,  I  is  the  identity 
matrix,  Ti  is  the  ith  column  of  si,  and  i  is  the  reference  sensor.  The  real 
scalar  s  is  the  signal-to-noise  ratio  (SNR)  to  be  used  for  filter  design. 

The  wavenumber  power  response  of  the  filter  F.  at  the  wave- 
number,  k  is  given  by 


where 


yt  =  C 


T  H  * 
V  F  F  V 
i  i  — 


j2TTk  •  xj  j2rtk  ■  x2 


j2rrk  •  x 


l\ 

/  and  Xj,  x2> 


,  xNare 


the  vector  distances  from  any  arbitrary  reference  point  and  the  N  sensors. 
The  sum  over  all  N  sensors  of  the  wavenumber  power  responses  of  the 
whitening  filters  designed  with  respect  to  each  sensor  is 


H  * 


T  *  1  _#  _T  T  *  T  *  T  * 

=  V  Z  V  A;  r.  r.  Z  V  =  V  Z  Z  V 

—  s2  — i  1  —  — 


Using  the  theorem  for  exact  inverse  of  A  +  B^C  and  sub- 
H  T 

stituting  si  =  A,  X  =  B^,  and  X  =  C)  ,  the  following  equation  is  obtained: 


Z  =  s 


(sI.XX11)'1  .  s|±.  -  1  =  I  -  -JUg- 

v  —  ;  [B  b2  +  sxhixJ  S  +  XHX 
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Then, 


T  *  T  *  # 
=  N-CV  X  (V  X  ) 


Note  that  the  constant  C  is  independent  of  the  wavenumber 

vector  k  and  thus  must  be  evaluated  only  once.  For  each  wavenumber  at 

which  the  power-response  sum  is  evaluated,  the  vector  dot  product  VTX* 

T 

must  be  evaluated  since  _V  is  a  function  of  wavenumber.  Multiplying  the 

T  * 

complex  number  V_  X  by  its  conjugate  and  by  C  and  then  subtracting  this 
real  number  from  N,  the  desired  response  sum  is  obtained.  Note  that  this 
response  sum  is  always  greater  than  0  (for  nonzero  s)  and  less  than  N.  There¬ 
fore,  its  reciprocal,  which  is  the  estimate  of  the  wavenumber  spectrum,  is 
always  positive  and  finite.  Note  also  the  lack  of  sensitivity  to  the  choice  of  s. 

B.  AVERAGE  OF  HIGH-RESOLUTION  WAVENUMBER  SPECTRA  COMPUTED 
WITH  DIFFERENT  REFERENCE  SENSORS 

If  the  sum  of  M  of  the  N  possible  individual  wavenumber  spectra 
is  desired,  this  is  given  by 


vTp  fhv* 
~  i  i  — 


M  sn 


which  is  not  as  easily  evaluated  as  was  the  reciprocal  of  the  sum  of  wave¬ 
number  responses.  However,  some  simplification  is  possible. 
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Let 


yT  F  FH  y*  _  yT  p  rT  *  * 

~  ._i  _J  -  _i  _i  —  yi  yi 


where 


=  r 

yi  y2  • '  *  yr 

1  =  vT 

xHx 

T  A  A 

L 

H 

J 

s  +  2  2£ 

Then,  the  sum  of  M  wavenumber  spectra  is  obtained  by  evaluating 


1 


y.  y. 

7i  7i 


Thus,  the  operations  required  to  evaluate  this  sum  at  each  wavenumber  are 
as  follows:  multiply  the  row  vector  V_  by  the  M  columns  of  the  constant 
matrix 


I  - 


-*H* 
s  +  XH  X 


corresponding  to  the  M  reference  sei.  *ors  selected;  multiply  the  resulting  M 
complex  numbers  by  their  conjugates;  and  take  the  reciprocal  of  each  re¬ 
sultant  real  number  and  sum. 

The  required  number  of  complex  multiplies  and  adds  isM(N+  1), 
which  is  well  above  the  N  +  1  needed  for  the  first  technique  but  well  below  the 
number  required  to  solve  for  each  filter  and  filter  response  individually. 


sclsno*  ssrvlcss  division 


The  second  technique,  although  requiring  vnore  computation, 
is  preferred.  When  summing  power  responses  over  all  N  multichannel  fil¬ 
ters  as  shown  in  the  first  technique,  the  spectral  window  effects  of  the  con¬ 
ventional  technique  are  regenerated.  Examination  of  the  equation  for  this 

sum  reveals  it  to  be  a  nonlinear  transformation  of  the  conventional  wave- 

T  #  T  * 

number  spectrum  given  by  3£  2£  V  .  ^  course,  summing  over  all  N 

filter  responses  is  not  necessary  and  a  direct  solution  similarly  derived  is 
available  for  the  sum  of  any  subset  of  N  filter  responses. 
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